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ABSTRACT 

We present a detailed model atmosphere analysis of two of the oldest stars 
known in the solar neighborhood, the high proper motion white dwarfs SDSS 
J110217.48+411315.4 (hereafter J1102) and WD 0346+246 (hereafter WD0346). We 
present trigonometric parallax observations of J1102, which places it at a distance of 
only 33.7 ± 2.0 pc. Based on the state of the art model atmospheres, optical, near-, 
mid-infrared photometry, and distances, we constrain the temperatures, atmospheric 
compositions, masses, and ages for both stars. J1102 is an 11 Gyr old (white dwarf 
plus main-sequence age), 0.62 Mq white dwarf with a pure H atmosphere and T c s = 
3830 K. WD0346 is an 11.5 Gyr old, 0.77 M Q white dwarf with a mixed H/He at- 
mosphere and T c g = 3650 K. Both stars display halo kinematics and their ages agree 
remarkably well with the ages of the nearest globular clusters, M4 and NGC 6397. 
J1102 and WD0346 are the closest examples of the oldest halo stars that we know of. 

Key words: stars: distances — stars: Population II — stars: atmospheres — stars: 
individual (SDSS J110217.48+411315.4, WD 0346+246) — white dwarfs 



1 INTRODUCTION 

The majority of the stars with an initial mass less than 
about 8 Mq end up as white dwarfs (WDs) . The short main- 
sequence lifetimes of intermediate-mass stars means that 
they spend almost all of their lives as cooling WDs, radiat- 
ing away their residual thermal energy slowly. To first order 
there is a simple rela tion between the age and luminosity of 
a WD l|Mestellll95"2r i. A typical 0.6 M Q pure H atmosphere 
WD c ools down to 3 8 00 K in ab out 10 Gyr dFontaine et al.l 
l200ll ). IWinget et all (|l987r i and iLiebert et al.l (|l988l ) were 
the first ones to use the oldest WDs in the solar neigh- 
borhood to constrain the age of the Galactic disk. Further 
studies based on nearby high proper motion WDs demon- 
strate that the oldes t disk WDs are about 8 + 1.5 Gyr old 
i|Leggett et alj Il998l : lHarris et all 120061 : iKilic et all l2010bl . 
and references therein). 

WD cosmochronology also constrains the ages of the 
oldes t halo WDs in the neares t globular clusters or in the 
field. lHansen et all (12004120071 ) use 100+ orbit Hubble Space 
Telescope (HST) observations of M4 and NGC 6397 and 
derive ages of 12.1 (with a 95% lower limit of 10.3 Gyr) 
and 11.5 Gyr, respectively. Another large HST program on 
47 Tucanae has recently been completed, but an age esti- 
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mate from the WD cooling sequence is not yet available 
l|Kalirai et all [2012). The extension of this method to halo 
WDs in the field is more problematic due to the unknown 
pop ulation membership of the high proper motion objects 
(see lOppenheimer et al.1 l2001bl : iBergeron et al.l 120051 ) . 

Fast moving halo WDs must exist in the solar neigh- 
borhood, though in relatively small numbers compared to 
the larger population of younger disk WDs. For example, 
the Besangon Galaxy model predicts 127 disk and 3 halo 
WDs p er square degree w ithin 1 kpc for a Galactic latitude 
of 45° (|Robin et al.ll2003l ). However, the halo WDs are pre- 
dicted to be fainter than V — 22 mag. Hence, deep, wide- 
field photometric and astrometric surveys are usually re- 
quired to identify a significant population of fiel d halo WDs. 
Such a survey based on the SDSS + USNO-B l|Munn et al.1 
2004) and the SDSS + Bok telescope proper motions iden- 
tified several spectroscop ically confirmed halo WD candi- 
dates (|Kilic et al. 2010a, bj). In addition, there are a dozen 
ultracool WDs that may be thick d isk or halo WDs (e.g., 
iGates et af1l2004l ; lHarris et al.ll2008l ). However, trigonomet- 
ric parallax measurements are required to confirm their halo 
membership. Large scale surveys like the Panoramic Survey 
Telescope & Rapid Response System (Pan-STARRS), Palo- 
mar Transient Factory, and the Large Synoptic Survey Tele- 
scope will be extremely use ful for the identific ation of large 
samples of faint halo WDs (|Tonrv et al.ll2012l ). 
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Figure 1. The trajectory of J1102 on the sky. Right panel: The 
same trajectory with the proper motion taken out. The tip of 
each arrow is the position from a single image, and the tail is the 
computed location based on the fitted trajectory including zero 
point, proper motion, and parallax. 



Even though the probability of finding a relatively 
nearby halo WD is small for any small patch of the sky, 
serendipitious discoveries of nearby halo WDs have hap- 
pened. WD0346 and J1 1 02 are the best examples of such dis- 
coveries. lHamblv et all l| 19971 ) identified WD0346 as a high 
proper motion WD (1.3" yr" 1 ) in the UK Schmidt Tele- 
scope photographic pla tes taken from 1987 to 1994. Sim- 
ilarly, lHall et all (|2008l l identified J 1102 as a high proper 
motion WD (1.75" yr -1 ) while searching for high-redshift 
quasar candidates in the SDSS spectroscopy data (the SDSS 
spectrum of J1102 shows a br oad emission feat ure at 5750- 
5900 A, which is an artifact). H all et al. | (|2008l ) explain the 
spectral energy distribution (SED) of J1102 with a 3830 K 
pure H atmosphere WD. However, without a parallax mea- 
surement, the mass, cooling age, and the kinematic mem- 
bership of this WD cannot be constrained accurately. 

Here we present parallax and mid-infrared photomet- 
ric observations of J1102, and a detailed model atmosphere 
analysis of both stars using parallax, optical, and infrared 
photometry. Our observations are discussed in Section 2, 
whereas our model atmosphere analysis and the nature of 
these objects are discussed in Section 3 and 4, respectively. 



For J1102, we measure a proper motion of fi = 
(-104.8 ± 0.9, -1741.8 ± 0.9) mas yr" 1 , and a relative par- 
allax and formal error of 28.8 ± 1.4 mas. Using the colors 
and magnitudes of the reference-frame stars, and increasing 
the uncertainty slightly to allow for unmodeled systemat- 
ics, we estimate the absolute parallax to be 29.6 ± 1.7 mas. 
Our proper motion measurement is relative to the chosen 
reference stars. Corrections to absolute proper motions are 
generally of order 10 mas yr" 1 . Our measurement is consis- 
tent with fj, = (—10 6.9, —1750.0) mas y r -1 in the absolute 
reference frame from iMunn et al.1 |2004l ). 

Figure 1 displays the trajectory of J 1102 on the sky, and 
the same trajectory with proper motion component taken 
out. This figure shows that our observations cover a large 
range of parallax factor for J1102, and the parallax is well 
co nstrained. We use t he full Bayesian formalism described 
in iThorstensenl (|2003f l to estimate the distance, including 
the Lutz-Kelker correction and prior information from the 
proper motion and very liberal limits on the luminosity. 
J1102 is only 33.7 ± 2.0 pc away from us. 

2.2 Mid-Infrared Photometry 

We obtained S pitzer InfraRed Array Camera (IRAC, 
iFazio et all 12004! ) 3.6, 4.5, 5.8, and 7.9 fim images of 
WD0346 and J1102 as part of the Cycle 3 program 30208 
and Cycle 5 program 474. For each object, we obtained 100 
second exposures for five dither positions in each filter. Our 
re duction procedures are similar to the procedures employed 
bv lKilic et all (|2009bh . Since our targets are relatively faint, 
we use the smallest aperture (two pixels) for which there 
are published aperture corrections. For J1102 we measure 
58.1±2.3, 40.4T3.0, 35.1T9.8, 23.3T13.9 ^Jy in channels 1, 
2, 3, and 4, respectively. Unfortunately, WD0346 is blended 
with a brighter source in the IRAC images, which prohibits 
us from performing aperture photometry o n it. None of 
these objects are detected in the WISE survey (|Wright et al.l 
l2010t ). 



2 OBSERVATIONS 
2.1 Parallax 

We obtained trigonometric parallax observations of J1102 
at the MDM 2.4m Hiltner telescope. We used a total of 106 
images taken on 13 observing runs between 2008 February 
and 2011 December. The instrumentation, observing pro- 
tocols, and reducti on techniques used w ere very similar to 
those descr ibed in Thorstcnscn ( 2003), iThorstensen et all 
(2008), and iKilic et al.l ( 2008aT ). For nearly all observations 
we used a 2048 2 SITe CCD detector ('Echelle') at the f7.5 
focus, but for the last two runs we used a 2048 2 STIS CCD 
('Nellie'), due to problems with the Echelle CCD. Echelle 
has 24 /tm pixels subtending 0.28", and Nellie has 21 fim 
pixels subtending 0.24". At each epoch we took many expo- 
sures in a 4-inch wide Kron-Cousins J-band filter, as near to 
the meridian as we could to minimize differential color re- 
fraction effects. The parallax reduction and analysis pipeline 
was unchanged from the previous work, and the change of 
detectors appeared to have no effect on the results. 



3 MODEL ATMOSPHERE ANALYSIS 

We use the state of the art WD model atmospheres 
to fit the available photometry for our targets. The 
model atmospheres includ e the Ly-a red wing opacity 
|Kowalski fc Saumonll2006T ) as well as non- ideal physics of 
dense helium that includes r efraction ([K owalski fc Saumonl 
l2004h . ionization equilibrium l|Kowalski et al .1120071) and the 
non-ideal dissociation equilibrium of H2 (|Kowalski|[2006al ). 
We convert the magnitudes into monochromatic fluxes us- 
ing the zero points derived from the Vega (STIS) spectrum 
integrated over the passband for each filter. We perform a 
two dimensional least squares fit in T e ff and the He/H ratio 
by minimizing the differences between the synthetic and the 
observed fluxes weighted by the observational errors. We use 
the parallax measurements to constrain the surface gravity. 



3.1 J1102+4113 

We u se the optical and near-infrared photometry (|Hall et al.l 
2008) plus our Spitzer photometry and distance measure- 
ment to model the J1102 SED. Figure 2 presents our best- 



11-12 Gyr Old WDs 3 




Wavelength (^m) 

Figure 2. The J1102 SED compared to the best-fit pure hydrogen 
atmosphere model spectrum (T c ff = 3830 K, solid line). Two 
other models with T off = 3600 K (dotted line) and 4000 K (dashed 
line) are shown for comparison. The optical spectrum is from the 
SDSS, with the artifact at 5750-5900 A blocked out. 

fit model compared to the observations. J1102 is best ex- 
plained by a r c ff = 3830 K, log g = 8.08, pure hydrogen 
atmosphere WD with M = 0.62M(T), which confirms the 
temperature assignment of lHall et al.1 (|2008l ). The cooling 
age of the WD is 10.015'^ Gyr. The error in the tempera- 
ture estimate is about 200 K. The best-fit model matches 
the infrared data fairly well; even though the H/He ratio 
is left as a free parameter in our model fits, the best-fit 
model always has a pure hydrogen atmosphere composition 
for this star. The fit to the optical portion of the SED is rea- 
sonable, but not perfect, possibly due to the problems with 
the Ly-a opacity calculations . The Ly-a opacity calculations 
from lRohrmann et al. I (|201ll ) do not result in any significant 
differences in our model calculations for this star. The op- 
tical portion of the SED favors a slightly cooler solution, 
whereas the infrared SED favors a hotter solution. There 
are known problems with the H2-H2 and H2-He collision- 
induced absorption (CIA. iFrommhold et al.ll201fj ) calcula- 
tions in dense WD atmospheres. Even though the strong 
CIA feature around 2/im is reproduced in our models, the 
predicted feature around has never been observed in 

cool WDs. Regardless of these problems, our best-fit model 
matches the overall SED fairly well. J1102 is clearly a very 
cool and old WD with a pure hydrogen atmosphere. 

3.2 WD 0346+246 

lHamblv et all l| 19991 ) measure a parallax of 36 ± 5 mas for 
WD0346, placing it at 28 ± 4 pc. The optical to near- 
infrared SED of WD0346 is kindly made available to us 
by B. Oppenheimer and is presented in Figure 3. Using 
rather crude mode l s for mixed H/He atmosphere WDs, 
Oopcnheimer et al l l|2001al ) find a best-fit solution of T e fr = 
3750 K and a He-dominated atmosphere with trace amounts 
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Figure 3. The WD0346 SED l|Qppenheimer et alJlioOlal l com- 
pared to the best-fit model with T e ff = 3650 K and He/H=0.43. 
Two other models with T e g = 3500 K (pure H atmosphere, dotted 
line) and T e ff = 3600 K (He/H=l, dashed line) are also shown. 



of hydrogen (logHe/H = 9.1). IBergeronl (|200ll ) points out 
that an extremely helium-rich composition is unlikely for 
such an old WD, because of accretion from the interstel- 
lar medium. In order to reproduce the short w avelength 
end o f the cool WD spectra, including WD0346, IBergeronl 
(2001) adds a pseudocontinuum opacity (due to the bound- 
free opacity associated with the so-called dissolved atomic 
levels of the hydrogen atom). The spectrum of WD0346 
is then fit fairly well with a T eS = 3 780 K, log g = 8.34 
and H e/H = 1.3 model. Mor e recently. iKowalski fc Saumonl 
(20061) and iKowalskH |2006b) identified the missing opacity 
in the blue as the red wing of the Ly a absorption, and mod- 
els including this opacity successfully reproduce the short- 
wave l ength fluxes of many coo l WDs (IKowalski fc Saumonl 
|2006j;lKilic et al.ll2008bl,l2009al lb1. l2010al ; lDurant et al. II2012I ; 
iGiammichele et alj|2012l ). 

In this paper we use the models with the aforemen- 
tioned Ly-a opacity to derive the atmospheric parameters 
of WD0346. The observed SED of WD0346 is best ex- 
plained by a mixed atmosphere model with T a g — 3650 K, 
logg = 8.3, and He/H = 0.43, which is fairly co n sisten t 
with the hydrogen-rich solution found by IBergeronl (|200lh . 
Figure 3 shows this model plus two other models with dif- 
ferent H/He compositions. Unlike J1102, the pure hydro- 
gen atmosphere model cannot reproduce the observed SED 
of WD0346 in the near-infrared. This is an indication that 
helium is present in the atmosphere and that through colli- 
sions with hydrogen species it enhances the CIA opacity. The 
best-fit mixed H/He atmosphere model matches the optical 
portion of the SED extremely well; the temperature of the 
star is constrained accurately by the optical data. However, 
the fit to the near-infrared data is not perfect. As in J1102, 
the models predict a broad absorption feature around 1-1.2 
/jm that is not observed. Regardless of these issues, WD0346 
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is clearly a very cool and old WD. Our best-fit model implies 
a mass of 0.77 M© and a cooling age of 11.2±? J Gyr. 



4 DISCUSSION 
4.1 Total Ages 

lHall et al] |2008f) classify J1102 as either a pure-hydrogen 
atmosphere WD at T e g = 3830 K (the most probable 
solution) or a mixed H/He atmosphere WD with H/He 
« 10~ 5 and T cB = 3500 K. Thanks to our Spitzer pho- 
tometry and parallax measurement, we are now able to 
constrain the atmospheric composition, mass, and age of 
this WD. Our detailed model atmosphere analysis based 
on optical/infrared photometry and distance measurements 
for J1102 demonstrate that it is a 3830 K, 0.62 M© star 
with a WD cooling age of 10.01?;* Gyr. The initial-to-final 
mass relation for WDs in dicate that the progenitor of J1102 
was a 1.8-2.2 M© star (ICatalan et all 120081 ; iKalirai et all 
120091 ; IWilliams et al] I2009F with a main-sequence lifetime 
of 0.6-1.1 Gvr |Marigo et al.ll200St ). Hence the total age of 
this object is 10.6-11.1 Gyr. Similary, WD0346 is a 3650 



K, 0.77 Mq star with a WD cooling age of 11.2±5'| Gyr. 
The progenitor star was a 3.1-3.3 M© main- sequence star 
with a main-sequence lifetime of 240-270 Myr (|Marigo et al.l 
2008). The uncertainty in main-sequence age is larger than 
the range given here. However, this uncertainty is insignif- 
icant compared to the total age of the WD, which is 11.5 
Gyr for WD0346. Both WD0346 and J 1 102 are significantly 
older than the coo lest disk WDs known IjLeggett et al.ll 19981 ; 
iKilic et aUkoiObft . 
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Figure 4. The Galactic orbits of J1102 and WD0346 for the past 
1 Gyr (solid lines) and the next 1 Gyr (dashed lines). The current 
position of each WD is marked with a dot. 



the plane by as much as 2.1 kpc and its closest approach 
to the Galactic center occurs at a distance of R ~ 2.2 kpc. 
Even though WD0346 is only 28 pc away from us right now, 
it will also travel as far away as 16 kpc within the next 1 
Gyr. WD0346 is most likely a halo WD. 



4.2 Halo Membership 

At a Galactic latitude of +64°, J1102 is 50 pc above the 
Galactic plane. Assuming zero radial velocity, it has U = 
68 ± 5, V = -259 ± 16, and W = 50 ±3 km s" 1 with r espect 
to the local standard of rest |Schonrich et al.l [2010). The 
unknown radial velocity mostly affects the W velocity. For 
example, a change in the radial velocity from —50 to +50 
km s _1 corresponds to a change in W from 5 to 95 km s -1 
and V from -261 to -256 km s" 1 . J1102 lags behind the 
Galactic disk. Figure 4 plots the Galactic orbit of J1102 
for the pa st and the next 1 G yr, in a static disk-halo-bulge 
potential l|Kenvon et al.| [2008V J 1102 goes above and below 
the plane by as much as 5 kpc and its closest approach to 
the Galactic center occurs at a distance of R « 150 pc. Even 
though J1102 is only 33.7 pc away from us right now, it will 
travel as far away as 16 kpc within the next 1 Gyr. J1102 is 
clearly a halo WD. 

At a Galactic latitude of -23°, WD0346 is 9 pc above 
the plane. Assuming zero radial velocity, it has U = —2 ± 
7, V = -148 ± 22, and W = -55 ± 9 km s" 1 with respect 
to the local standard of reslQ. The unknown radial velocity 
mostly affects the U velocity. For example, a change in the 
radial velocity from —50 to +50 km s _1 corresponds to a 
change in U from +43 to -47 km s" 1 . WD0346 also lags 
behind the Galactic disk. WD0346 goes above and below 

1 The UVW vel ocities for WD0346 were miscalculated by 
Kilic et al.l j2010al ) due to a sign error in its proper motion. 



4.3 The Oldest Stars in the Solar Neighborhood 

We n ow know of star s that are almost as old as the Uni- 
verse |Hill et al.l l2002). However, old stars with reliable age 
measurements are rare; radioctive elements like thorium and 
uranium are detected only in a fe w cases, and even t hen the 
age uncertainties are ~ 2-3 Gyr (|Frebel et al.ll2007l ). Com- 
pared to main-sequence and giant stars, WDs provide an 
advantage in age measurements, thanks to the relatively 
short main-sequence lifetimes of typical 0.6 M© WDs and 
the w ell understood physics of the cooling WDs. ISion et all 
(2009) pres ent a volume limited samp le of WDs within 20 pc 
of the Sun. iGiammichele et al.l (|2012l ) show that the coolest 
WDs in that sample are about 4200 K; significantly hotter 
than the two WDs discussed in this paper. 

With estimated ages between 11 and 12 Gyr, J1102 
and WD0346 are currently the oldest stars known in the 
solar neighborhood. The theoretical uncertainties due to 
the unknown core composition, helium layer mass, crystal- 
lization, and phase sepa ration are ~1 Gyr for these ages 
(| Montgomery et al.lll999l ). The age estimates for these two 
WDs are remarkably similar to the oldest WDs found in 
the nearest globular clusters M4 and NGC 6397. J 1102 and 
WD0346 are about 100 times closer than these clusters. 
They provide an unprecedented opportunity to understand 
the model uncertainties in cool WD atmospheres and put the 
globular cluster ages on a more secure footing. Our analysis 
shows that the current WD atmosphere models including the 
Ly-a opacity in the blue and the CIA in the infrared provide 
a reasonable match to the observations, though problems 
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most likely associated with the CIA calculations remain. In- 
frared data are essential for constraining the atmospheric 
composition of cool WDs (see Fig. 2 and 3). Hence, the 
temperature and age estimates based on two filter HST op- 
tical data are likely uncertain by a few hundred degrees and 
several hundred Myr, respectively. 



5 CONCLUSIONS 

We present trigonotometric parallax and Spitzer mid- 
infrared photometric observations of the halo WD candi- 
date J1102. We use these new data along with optical and 
near-infrared data on J 1102 to constrain the temperature, 
mass, and age of this WD. J1102 is an 11 Gyr old WD 
33.7 pc away from the Sun. We also revisit the model at- 
mosphere analysis of another halo WD candidate, WD0346, 
using improved model atmosphere calculations. WD0346 is 
an 11.5 Gyr old WD 28 pc away from the Sun. WD0346 
and J 1102 are currently the oldest stars known in the solar 
neighborhood. Both stars display halo kinematics; they are 
just visiting our neighborhood at the moment. 

WD0346 and J1102 are the best examples of serendip- 
itious discoveries of nearby halo WDs. Current and future 
deep, wide-field surveys ought to find many more old halo 
WDs. Such WDs remain to be discovered in Pan-STARRS, 
Palomar Transient Factory, the Large Synoptic Survey Tele- 
scope, and GAIA data. These discoveries will provide inde- 
pendent age measurements and constrain the age and age 
range of the Galactic halo. 
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